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1. Introduction 

The CGC picture [1] of high energy QCD has been succes- 
fully applied to describe deep inelastic scattering, single particle 
production and particle correlations in forward rapidity proton- 
proton and proton nucleus collisions. It also forms a basis for 
understanding the initial stage of deconfined QCD matter cre¬ 
ated in ultrarelativistic nucleus-nucleus collisions. An essential 
ingredient in these calculations are the JIMWLK equation and its 
mean field limit — the Balitsky-Kovchegov (BK) equation [2, 3]. 

They describe the energy, or equivalently Bjorken-.r, dependence 
of correlators of Wilson lines in the target color field. Cross sec¬ 
tions for differenf processes are fhen expressed in ferms of fhese 
Wilson line correlafors. 

Recenfly firsf sfeps have been faken fo develop fhis phenomenological picfure fo nexf-fo- 
leading order (NLO) accuracy, bofh for fhe evolution equafions fhemselves [4, 5, 6] and for specific 
scattering processes [7, 8, 9, 10]. There remain bofh concepfual and numerical challenges in carry¬ 
ing ouf fhe full CGC phenomenology program af NLO accuracy. This falk discusses in more defail 
one of fhese aspecfs, by presenting fhe resulf [11] from a direcf “brufe force” numerical solufion of 
fhe NLO BK equation as if is wriffen down in Ref. [4]. 



Figure 1: Coordinates for the 
eikonal Wilson lines and their 
separations 


2. The equation 


In the large-Ac limit (but assuming Nf ~ Ac) we can write the evolution equation derived in 
Ref. [4] for the dipole as: 

dyS{r) = ® [S{X)S{Y)-S{r)] + ® S(F)[S(A') -5(X)] 

+ ^^2 ® [S{X)S{z-my')-S{X)S{Y)], with S(r) = 1 {TrU^{x)U{y)) , (2.1) 

Here U {x) is a fundamental representation Wilson line describing the propagation of a high energy 
probe through the dense color field of fhe fargef. The Wilson lines are needed af coordinafes 
x,y,z,z' in fhe two-dimensional transverse plane, with the six distances between them denoted as 
r,X,X',Y,Y' and z — z' as shown in Fig. 1 . The convolutions (g) denote integrations over z (in ) or 
z and z' (the other terms). We replace the terms explicitly proportional to the j8-function coefficient 
in the first kernel Ki by the “Balitsky” running coupling prescription [12]. With this replacement 
the explicit expression for the first kernel is 
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The first term inside first the square bracket in Ki is the leading order fixed coupling kernel. The 
whole firsf square brackef in Ki corresponds fo fhe “Balifsky” running coupling prescripfion oflen 
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used in LO phenomenology [13, 14], with the second square bracket being a pure NLO term. The 
coupling constant in front of the purely NLO kernels K 2 and Kf is taken to depend on the parent 
dipole size r and the explicit expressions for the kernels are 


2 lX^Y'^+X'^Y^-4r^{z-z')^ 

{z-z'Y^[ {z-z'f{X^Y'^-X'^Y^) 


^4 ^2-1 

^ X^Y'^{X^Y'^-X'^Y^)^ X^Y'^{z-z'Y\ ^X'2f2 

^ {z-z'Y {z-z'Y{X^Y'^-X'^Y^) X'2t2- 


(2.3) 


(2.4) 


In addition to the j3-function terms, two kinds of logarithms appear in the kernels. The ones 
in K 2 and Kf depend on conformal ratios of four distances, and vanish in the limit r —)■ 0. The 
first kernel Ki , on the other hand, has a nonconformal double logarithm that diverges in the limit 
r —)• 0. Although this is an integrable singularity in the z-integral, it nevertheless has a problematic 
effect on the evolution equation, as we will show in the following. As an initial condition we use a 
parametrization 


N{r) = 1 —S{r) = 1 — exp 





(2.5) 


with two tunable parameters: 


• The ratio 2 so/''^qcd essentially determines value of ots and controls the overall relative im¬ 
portance of the NLO corrections. 

• The anomalous dimension 7 controls the shape of the initial condition. Leading order fits to 
HERA data using the parametrization (2.5) prefer a value 7 > 1 which then becomes 7 ~ 0.8 
during the BK evolution. Since in an NLO fit also the impact factor relating the cross section 
and the dipole amplitude N{r) should be different from the LO one, without performing the 
full fit it is not a priori obvious what would be a value favored by experimental data. 


3. Properties of the solution 

Figure 2 (left) shows the logarithmic evolution speed dyN{r)/N{r) for the MV model initial 
condition 7 = 1 . Firstly it is obvious that for small values of 2s/AQCD^ i-S- effectively large cou¬ 
plings, the the evolution speed is negative at all values if r. This means that the NFO corrections are 
large and negative, and the scattering amplitude actually decreases with energy; a very nonintuitive 
result. For smaller typical values of tts the region around the “front” r ~ l/2s behaves in a rea¬ 
sonable way, but for very small dipoles the evolution speed still diverges as dyN/N ~ In r. Figure 2 
(right) shows the contributions from different parts of the equation. The FO term gives a positive 
dyN/N that approaches a constant at r —)• 0. The divergence for small r is due to the nonconformal 
double logarithm, while the other NFO corrections yield a contribution that is negative, but smaller 
in magnitude than the leading order result. 

While having dyN{r) < 0 or V(r) < 0 is not very physical, it is not in itself a mathematical 
contradiction. A divergent dyN /V ~ Inr in the limit r —)• 0 is, however, a signal of an instability in 
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Figure 2: Left; Evolution speed dyN{r)/N{r) at the initial condition y = 0 for the MV model initial condi¬ 
tion, 7=1 Right; Contribution of different parts of the equation (2.1) to the evolution speed dyN{r)/N{r). 
Also the dipole amplitude N{r) is shown with a thick grey line. 

the system. One way to see this is the following simple argument. Let us consider a small but finite 
interval in rapidity, dy, as in a numerical solution of the evolution equation. A diverging dyN/N 
for r —1 0 means that there is a small but finite r below which N becomes negative already in this 
one step in rapidity. This immediately makes the equation unstable, since the convergence of the 
z-integral of the leading order equation requires N{r) —0 for r —0. A finite N{r = 0) is also 
inconsistent with the definition A(;c — y) = 1 — ^TrU\x)U(y) in terms of Wilson lines. To avert 
this problem in the numerics, we always enforce A(r) > 0 by hand. 


4. The anomalous dimension 


<5.s,o/AqcD — 19 Qs.o/Aqcd — 19 Q,s.oMqcd — 19 



Figure 3: Anomalous dimension /(r) = dlnV(r)/dlnr^ aty = 1, 5, 30. 


Some more insight into the behavior of the equation can be obtained by following the devel¬ 
opment of the r-dependent anomalous dimension that we define here as 


7(r) 


dlnA(r) 

dlnr^ 


(4.1) 


In terms of the anomalous dimension we can discuss the behavior of the equation at small r by 
parametrizing the amplitude as N{r) ~ (Qsr)^^. For an appropriate 7 the leading order equation 
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maintains this form, with ~ and dyN/N ^ 2yX > 0. If, due to large NLO corrections, one 
has dyN jN —)> —c < 0 for r —)> 0, this means that the solution behaves as N ~ mathematically 
this is not problematic, but physically it is unnatural to have the scattering amplitude decrease 
with energy. If, however, dyN/N ~ clnr, we can parametrize N{r) ~ with /(y) ~ y. 

In other words, for a diverging evolution speed the functional form of the amplitude as a function 
of r gets steeper with the evolution. If one enforces N > 0 for small r and N < 1 for r —)■ oo, the 
amplitude eventually develops into a discontinuous form N(r) ~ 6(r — l/2s)- The evolution of 
7 (r) for different initial conditions is shown in Fig. 3, where this behavior can be clearly seen. The 
unstable nature of the equation shows up as an increase in y/r) for small r. This increase starts 
immediately for the initial condition with 7=1, more slowly for the initial condition 7 = 0.8 and 
is not noticeable within the rapidity range studied here for y = 0.6. 


5. The composite conformal dipole 


In Ref. [15] Balitsky and Chirilli note that the nonconformal double logarithmic term appears 
also in the equation for the conformal N = 4 Super Yang-Mills theory. It is therefore interpreted 
as an artefact of a cutoff that breaks conformal invariance. The authors then propose to correct 
for this effect order by order in perturbation theory by introducing a “composite conformal dipole” 
operator defined as 


= S{r) - 


«sfVc 

4k'^ 


d/z 


In 


ar 


X2y2 


-[SiX)SiY)-S{r)], 


(5.1) 


Qsfl/ Aqcd = 19 


with a a dimensionful constant that drops out 
of the final equafion. Rewrifing fhe equafion 
(2.1) in terms of fhe conformal dipole (5.1) 
removes fhe nonconformal double logarifm, 
buf infroduces an addifional term in fhe ofher 
kernel K 2 fhaf behaves as Inr^ for small r. 

We have checked numerically fhaf fhe quali- 
fafive behavior of fhe conformal dipole equa¬ 
fion remains fhe same as fhe original one. As 
shown in Fig. 4, if is now fhis new Inr^ term 
which is responsible for fhe leading behavior 
af small r. 

In conclusion, we have performed fhe 
firsl numerical solufion of fhe full NLO BK 
equation direcfly in coordinate space. The 
NLO corrections are negafive, implying a 

slower energy dependence of cross secfions fhan wifh fhe LO equation. This generically leads 
fo a heifer agreemenl wifh experimenlal dala. The equafion has, however, a double logarilmic lerm 
fhaf causes a problemalic behavior for small dipoles, i.e. high Q^. If seems evidenl fhaf Ihese large 
logaiilhms will need fo be resummed in order for fhe equafion lo be useful in praclical phenomeno¬ 
logical work. 



Figure 4: Contribution of different terms to the evo¬ 
lution speed dyN/N of the conformal dipole. 
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